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Abstract. Hard scattered partons are predicted to be well calibrated probes of the 
hot and dense medium produced in heavy ion collisions. Interactions of these partons 
with the medium will result in modifications of internal jet structure in Au+Au events 
compared to that observed in the p+p/d+Au reference. Full jet reconstruction is 
a promising tool to measure these effects without the significant biases present in 
measurements with high-px hadrons. 

One of the most significant challenges for jet reconstruction in the heavy ion 
environment comes from the correct characterization of the background fluctuations. 
The jet momentum irresolution due to background fluctuations has to be understood 
in order to recover the correct jet spectrum. Recent progress in jet reconstruction 
methodology is discussed, as well as recent measurements from p+p, d+Au and Au+Au 
collisions at ^snn = 200 GeV. 



1. Introduction 

Jets are remnants of hard-scattered partons, which are the fundamental objects of 
perturbative QCD. At Relativistic Heavy Ion Collider (RHIC), they can be used as a 
probe of the hot and dense matter created in heavy ion collisions. Interaction and energy 
loss of energetic partons in the medium lead to jet quenching in heavy ion collisions. 
Until recently, jet quenching was studied indirectly using single particle spectra and 
di-hadron correlations [1]. These measurements are however limited in the sensitivity 
to probe partonic energy loss mechanisms due to biases toward hard fragmentation and 
small energy loss |2]. 

Developments in theory (for example [31 H]) and experiment (detector upgrades, 
increased RHIC luminosity) finally enabled full jet reconstruction in heavy ion 
collisions |5]. Full jet reconstruction reduces the biases of indirect measurements and 
enables access to qualitatively new observables such as energy flow and fragmentation 
functions. As a baseline measurement for heavy ion jet studies, p+p collisions at 
the same energy are used. To isolate initial state effects from medium modiflcation, 
measurements in d+Au are essential. 

We present current jet analyses at STAR, starting with recent results on initial 
state effects (d+Au). Status of jet spectra analysis in Au+Au follows, including studies 
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of background fluctuations and their effect on the measurement of jet spectra. Finally 
we discuss recent results on correlations in Au+Au triggered by fully reconstructed jets: 
di-jet analysis and jet-hadron correlations. 

2. Jet reconstruction 

The present analysis is based on ^snn = 200 GeV data from the STAR experiment, 
recorded during 2006-2008. The Barrel Electromagnetic Calorimeter (BEMC) detector 
is used to measure the neutral component of jets, and the Time Projection Chamber 
(TPC) detector is used to measure the charged particle component of jets. In the 
case of a TPC track pointing to a BEMC tower, its momentum is subtracted from the 
tower energy to avoid double counting (electrons, MIP and possible hadron showers in 
the BEMC). Pseudorapidity acceptance for jets is \ri\ < 0.6 in p+p and Au+Au and 
\ri\ < 0.55 in the case of d+Au colhsions. 

Recombination jet algorithms and anti— /ct, part of the Fast Jet package [3], are 
used for jet reconstruction. To subtract the background, a method based on active jet 
areas [4] is applied event-wise: p^'^'^ = p^^^'^^'^^^'^ — p- A, with p estimating the background 
density per event and A being the jet active area. 

An important aspect of underlying event background are its fluctuations. We 
discuss data-driven methods used to correct the jet observables for these fluctuations. 

3. Initial state: d-|-Au 

This analysis is based on minimum bias triggered ^snn = 200 GeV data from the 
STAR experiment, recorded during RHIC run 8 (2007-2008). The Beam Beam Counter 
detector, located in the Au nucleus fragmentation region, was used to select the 20% 
highest multiplicity events in d+Au collisions. lOM events after event cuts were used for 
jet finding (anti— /ct algorithm) with a resolution parameter R = 0.4 andpx > 0.2 GeV/c 
cut was applied to tracks and towers. 

PYTHIA 6.410 and GEANT detector simulations (adjusted to match the realistic 
TPC tracking efficiency in d+Au run 8 running) were used for jet corrections to hadron 
level. Embedding into real d+Au events at level of reconstructed tracks and towers was 
used to correct for background fiuctuations. A bin-by-bin correction was applied to the 
jet spectrum [_6J. 

To compare the per event jet yield in d+Au to jet cross section measurements in 
p+p collisions, MC Glauber studies were utilized: (Nun) = 14.6 + 1.7 for 0-20% highest 
multiplicity d+Au colhsions and crinci,pp = 42 mb. These factors were used to scale the 
p+p jet cross section measured previously by the STAR collaboration |7] using a Mid 
Point Cone (MPC) jet algorithm with a cone radius of R = 0.4. The resulting d+Au jet 
Pt spectrum is shown in Figure [1] together with the scaled p+p jet spectrum. Within 
the systematic uncertainties, the d+Au jet spectrum scales with (A^bin)- 

The leading systematic uncertainty is the Jet Energy Scale (JES) that is driven 
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by imprecise knowledge of TPC tracking efficiency for tracks in jets with realistic run 
8 d+Au detector backgrounds. This will be considerably improved by embedding jets 
in raw d+Au data. With better handle on JES and by measuring jet spectrum in run 
8 p+p collisions and in peripheral d+Au collisions, we'll be able to construct Raa and 
Rcp for jets, respectively. 



4. Inclusive jet spectra and background fluctuations in Au+Au 

Preliminary results on jet spectrum in Au+Au collisions at ^snn = 200 GeV 
were reported in [S]. In this analysis, the background fluctuations were estimated by 
generating PYTHIA jets and embedding them into real central Au+Au events. The 
resulting spectrum distortion was parametrized by Gaussian, for R = 0.4 the width 
is 0" = 6.8 GeV with systematic uncertainty ±1 GeV. This parametrization was then 
used for a regularized matrix inversion to unfold the measured jet spectrum. Resulting 
Raa is shown in Figure El The systematic uncertainties prevent us from precisely 
quantifying the suppression for R = 0.4 jets (their reduction will be the subject of the 
next paragraphs). It's clear however, that these jets are less suppressed than jets with 
R = 0.2 and charged hadrons (-Raa ~ 0.2). This is consistent with a picture of jet 
proffie broadening from R = 0.2 to R = 0.4 in central Au+Au collisions with respect to 
p+p, which is illustrated by spectra ratios in Figure [31 
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Figure 1. Jet px spectrum from 
d+Au collisions 6 compared to 
(-^bin) scaled p+p spectrum 
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Figure 2. Raa of jets in 
central Au+Au collisions for 
and anti— fcx algorithms 
and i? = 0.2, 0.4 [8]. 



Precise characterization of underlying event background fluctuations is essential 
to reduce systematic uncertainties in jet measurements. These are hence a subject 
of intense study, both theoreticUy [9] and experimentally. We summarize here recent 
results of STAR studies of background fluctuations [TOj . 

To quantify the background fluctuation, a method is used where a probe "jet" 
(single particle, PYTHIA jet, QPYTHIA jet) is embedded into real central Au+Au 
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events. This allows to extract the shape of the fluctuations (they are known to be non- 
Gaussian) and also to check if they are independent of the fragmentation pattern of 
the probe (this is essential as jet fragmentation is modifled in Au+Au collisions). We 
embed an object of known = p^^^'^ and apply jet reconstruction on the hybrid event 
(anti— /ct algorithm with R = 0.4). We match a reconstructed jet containing more than 
50% of probe to the probe jet and quantify the response of the hybrid system to the 
embedded jet via: 

Sp^ = - p ■ A''^™ - plp'^^*'^ (1) 

where A^'^'^° is the area of the matched reconstructed jet and p is determined prior 
to the embedding step. This deflnition is identical to Eq. (1) in [9J. The normalized 
distribution of SpT is the probability distribution to flnd jet energy (after event-wise 
background correction) p^^^ = pt^'^^ + Sp^^. If there were no background fluctuations, 
5pT would be a delta function at zero. For very low px probes, areas of anti— /ct jets get 
very small, so a cut A^'^'^" > 0.4 was apphed. With this cut, 5pT distribution turns out 
to be largely independent of p^^^"^ [10] . 

We have investigated dependence of (5pT on jet fragmentation pattern. Figure H] 
shows the overlay of multiple 5pi distributions for single particle jets and for jets with 
both low and high px generated by PYTHIA and Q-PYTHIA (g = 5 GeVVfm). In 
order to compare their shapes directly, the distributions were aligned horizontally by 
fitting a Gaussian function to 5p-£ < and aligning the centroids by shifting relative to 
one reference distribution. The shifts are shown in the insert and are typically smaller 
in magnitude than 1 GeV. The overlay shows that the 6pi distribution is to a large 
extent universal, within a factor ~ 2 at 5pi = 30 GeV, especially in region 6pT > 
which drives the smearing of the inclusive jet spectrum. Further quantification of this 
observation and its application to deconvolution of the measured inclusive jet spectrum 
in central Au+Au collisions is in progress. 

5. Jet triggered correlations 

A highly biased jet population was used as trigger in di-jet and jet-hadron correlations. 
Trigger jets are required to contain a BEMC tower with Et > 5.4 GeV to achieve a 
longer in- medium pathlength on the recoil (away side). To limit the effect of background 
fluctuations, trigger jets are reconstructed using only TPC tracks and BEMC towers 
with Pt > 2 GeV/c. A 2 GeV systematic uncertainty on trigger jet energy was used to 
account for any remaining effect. 

Recoil jet px spectrum was measured in p+p and Au+Au (0-20% most central) 
collisions [11]. A Gaussian model of background fluctuations was used to unfold the 
Au+Au spectrum with systematic uncertainty ±1 GeV. Figure [5] shows a signiflcant 
suppression of recoil jet pt spectrum in Au+Au compared to p+p for R = 0.4, which 
suggests jet broadening beyond R = 0.4. Considering also the observation from inclusive 
jet analysis (suggestive of jet broadening from R = 0.2 to i? = 0.4) there appears 
to be a smooth jet broadening trend. Note that the recoil jet px spectrum is much 
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Figure 3. Ratio of i? = Figure 4. Quantifying tlie back- 

2/R = 4 jet pt spectra ground fluctuations and their depen- 

in p+p and Au+Au colli- dence on probe px and fragmenta- 
sions 0. ElO]- 



flatter (harder) than the inclusive one, therefore impact of uncertainties on background 
fluctuations is much reduced. 

Full jet reconstruction is not feasible for R > 0.4 due to large background 
fluctuations. To investigate the jet broadening on the away side of the trigger jet, jet- 
hadron (JH) azimuthal correlations between trigger jet and charged hadrons (detected 
by the TPC detector) are measured [T2]. The raw azimuthal correlation is parametrized 
via two Gaussian peaks (near and away side jet) and V2 modulated background (with 
fixed V2 values). The uncertainties in the (a priori unknown) jet V2 value were chosen to 
cover the extreme cases of no V2 and 50% higher than f2{2} at px = 6 GeV/c (default 
is V2{2} at Pt = 6 GeV/c). The associated track V2 values and uncertainties follow the 
analysis in [13]. Due to ambiguities of ZYAM for broad jet structures, the background 
level was determined by the fit. For comparison ZYAM was applied (as expected for 
broad structures it leads to an underestimation of the correlated away-side yields for 
lower associated px)- 

Figure |6] shows the awayside Gaussian width of JH in p+p and — 20% most 
central Au+Au collisions. There is a significant broadening (Au+Au w.r.t. p+p) for 
p^ssoc ^ 2 GeV/c, while no broadening at higher observed. Jaa, the ratio of 

per-trigger associated yields, is plotted in Figure [71 There is a significant suppression 
of high Pt particles on the away side accompanied by an enhancement at low p^^°'^. In 
order to quantify the energy redistribution on the away side, it's better to instead of 
Jaa use Daa- 

DAAiprn = pT"" ■ (YAAiprn - (2) 

where Xaa,pp are per-trigger associated yields in AA,pp. Away side Daa for JH is 
shown in Figure [8l In fact, the energy "lost" at high pt is approximately compensated 
by low Pt enhancement [11]: jet quenching in action. 
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Given the observed broadening for p^^**"^ < 3 GeV / c, no broadening for high p^^"'^ 
and /aa shape independent of p^'^"'^ at high p^'"^"'^, one can speculate that the original 
parton loses energy by emission of soft radiation (and therefore the original jet direction 
changes little: no broadening is observed at high p^^""^). These soft fragments traverse 
the medium, receive transverse kicks and therefore appear at large angles with respect 
to the original parton direction. The energy loss is followed by a possibly vacuum-like 
fragmentation of a parton with reduced energy. 
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Figure 8. Away side Daa in JH 
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6. Summary 

We have presented STAR results on full jet reconstruction in d+Au collisions. Within 
current systematics there appears to be binary collision scaling compared to p+p 
collisions, but the final measurements (jet Raa and Rcp) with reduced systematic 
uncertainties are yet to be completed. The study of background fluctuations 6pT in 
Au+Au collisions suggests its independence of p^™tied ^j^g probe jet fragmentation 
pattern. The shape of 6pT will be used to unfold the measured jet px spectrum to obtain 
the final result with decrased systematic uncertainties. The hints of jet broadening 
obtained first in inclusive jet p^ spectrum and di-jet correlations were further studied 
using jet-hadron correlations, where significant broadening and enhancement at low 
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p^soc ^g^g accompanied by suppression (and no broadening) at high p^^"*^. This is 
consistent with the picture of pQCD-hke energy loss to soft fragments followed by 
(possibly vacuum-like) fragmentation of a parton with reduced energy. 
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